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FOREWORD 


This  report  presents  a  mathematical  derivation  of  the  relativistic  effects  in  satellite  tracking. 
Although  the  results  are  not  new,  the  presentation  of  the  derivations  in  this  systematic,  self-contained 
way  can  help  the  user  easily  make  modifications  to  suit  his  requirements. 

This  report  was  reviewed  by  Mr.  Robert  W.  Hill,  Head  of  the  Space  Flight  Sciences  Branch,  Space  and 
Surface  Systems  Division. 
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INTRODUCTION 


Since  the  advent  of  the  NAVSTAR  GPS  (Global  Positioning  System)  era,  it  has  become  more 
important  to  correct  for  relativistic  effects  in  satellite  orbit  determination.  This  is  primarily  so  because 
range  to  the  GPS  satellites  is  now  being  determined  by  timing  the  propagation  of  a  pulse  from  the 
satellites  to  the  receivers.  Prior  to  GPS,  Doppler  data  were  used  to  determine  the  orbits  from  range 
difference  or  range  rate. 

When  using  Doppler  data,  the  relativistic  effects  can  be  absorbed  into  a  frequency  bias  that  is 
formally  eliminated.  No  explicit  mathematical  equation  needs  to  be  incorporated  into  the  computer 
program.  This  assumes,  of  course,  that  the  periodic  relativistic  effect  due  to  orbital  eccentricity  is 
negligible. 

When  using  range  data,  however,  the  situation  becomes  more  complicated.  A  pulse  emitted  from 
the  satellite  carries  information  telling  the  proper  time  of  the  transmitter  (satellite)  clock  at  the  time  of 
emission,  rTg .  When  the  pulse  is  received,  the  proper  time  of  reception  is  noted  by  the  receiver  (station) 
clock,  rRr.  The  relativistically  uncorrected  range  then  is 


P u  =  c^Rr  -  TTe) 

where  c  is  the  vacuum  speed  of  light.  We  say  this  range  is  relativistically  uncorrected  because  two 
clocks  at  different  gravitational  potentials  do  not  run  at  the  same  rate:  for  two  clocks  at  rest,  the  one 
in  the  stronger  gravitational  field  runs  slower. 

Because  the  clocks  run  at  different  rates,  even  if  they  are  perfect  clocks  and  are  perfectly  synchro¬ 
nized  initially,  they  will  grow  more  and  more  out  of  synchronization  and  rre  ^  rRe  where  rRg  is  the 
proper  time  of  emission  as  read  on  the  receiver  clock.  The  relativistically  corrected  range  then  is 


P  =  c^Rr  ~  tr). 


Although  the  equations  for  the  relativistic  corrections  are  well  known  and  have  been  published  in 
several  papers,1,2  the  purpose  of  this  paper  is  to  derive  them  in  a  systematic  way.  The  only  part  that  is 
not  derived  is  the  Schwarzschild  solution  to  Einstein’s  field  equations.  All  other  derivations  are  in¬ 
cluded  in  this  paper  in  addition  to  the  Doppler  effect. 
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PRELIMINARY  DISCUSSION 


From  the  field  equations  of  general  relativity,  the  squared  line  element  in  a  centrally  symmetric 
gravitational  field  is 


ds2  =  - 


/  2GM\ 1 

1 


dr 2  -  r2  dO 2  -  r2  sin2  0  dip2  +  c2  [  1  - 


2GM 

c2r 


\dt 2 


(1) 


where  G  is  the  gravitational  constant,  M  the  mass  of  the  central  body,  and  c  the  velocity  of  light  in 
vacuum.  This  exact  solution  to  the  field  equations  was  obtained  by  Schwarzchild.3,4 

To  first  order  in  1/c2,  Equation  1  can  be  written  as 

ds 2  -  (  1  +  \dr 2  -  r 2  dd2  -  r1  sin2  dd<j>2  +  c2 


=  c2  dt 2 


=s  -  ^1  +  dr 2  -r2  dd 2  -  r1  sin2  dd<p2  +  c2  ^1  -  — y— ^dt2 

i/iy  .4/ay .  4  sin. eM  -  +  (,  -  ^ 

c^Xdt )  c‘\dtj  c2  \dt)  c4r\dtj  ^  c2r  J 


(2) 


The  first  three  terms  in  the  bracket  are  -  expressed  in  spherical  coordinates  and  since 
u2  2  GM  (  dr\2 

^  (  —  I  ,  the  fourth  term  can  be  neglected  and  Equation  2  can  be  written  as 


c*r  \dtj 


ds  =  cdt 


2GM  v2 


c2r 


(3) 


The  Newtonian  potential  for  a  spherically  symmetric  field  due  to  a  sphere  or  point  mass  is 

n 

4»(r)  = - (ji  =  GM)  so  in  terms  of  the  Newtonian  potential,  Equation  3  can  be  written 

r 


ds  =  cdt 


1  + 


2  $(/■)  u2 


(4) 


Proper  time  dr  is  the  time  interval  measured  by  an  observer  at  rest  in  any  system  and  is  defined  by 


ds 

dr  = —  . 

c 


An  observer  at  rest  in  a  system  that  is  moving  with  a  velocity,  v,  in  a  potential  4>(r)  then  will  measure 
a  proper  time  interval  given  by  Equation  4 


r  .  2<P(r)  u2> 
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The  quantity  dt  is  called  the  “coordinate  time  interval.”  For  two  observers  with  different  velocities 
and  at  different  potentials,  t/r,  =£dr2  whereas  dtx  =  dt2\  i.e.,  dt t  =  dt2  as  objects  observed,  not  as 
observers. 

If  we  imagine  all  points  in  space  containing  clocks  which  by  some  mechanism  are  kept  in  sychroni- 
zation,  the  clocks  will  always  be  in  agreement  and  show  the  coordinate  time.  This  coordinate  time  will 
be  the  same  as  the  proper  time  for  an  observer  at  rest  at  infinity  since  from  Equation  5  as  v~*  0 
and  0,  dr~ *  dt. 

The  distinction  between  proper  time  and  coordinate  time  is  very  important.  To  paraphrase 
Eddington,4  proper  time  is  the  time  of  consciousness  and  coordinate  time  is  the  time  in  physical  and 
astronomical  reckoning. 


RELATIVISTIC  CORRECTION  TO  RANGE 


SATELLITE  TO  EARTH  STATION 


We  can  write  Equation  5  for  a  transmitter  (any  satellite,  subscript  T)  and  for  a  receiver  (earth 
fixed  station  or  any  other  satellite,  subscript  R)  and  we  have 


drT  -  J^l  + 

Jtr  =  [> 


2<I >(rT) 

c 2 

¥  * 

(6) 

2<t>(rR ) 

IT  * 

(7) 

c2 

C2J 

Note  that  there  is  no  subscript  on  the  coordinate  time  interval,  dt,  as  discussed  in  the  previous  section. 
The  reason  for  no  subscripts  is  that  we  are  assuming  that  the  transmission  time  has  already  been  taken 
into  account.  We  are  only  comparing  clock  rates. 

When  Equation  6  is  divided  by  Equation  7 


which  gives  to  order 


(8) 


■J 

4 


4 


1 


1 


3 
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Potential  energy  is 


V  =  m<t>(rT) 


and  the  nonrelativistic  kinetic  energy  (valid  to  our  order  of  approximation)  is 


K  =  -  /huL 
2  T 


Also,  from  classical  mechanics  the  satellite’s  total  energy  is 


E  =  K  +  V 


pm 

57 


where  n  =  G/V/,  m  is  the  satellite  mass,  and  aT  is  the  semimajor  axis.  From  above  then 


um 

K  =  E  -  V  =  -  - - V 

2ar 


and  the  first  bracket  in  F.quation  8  becomes 


[* 


„2 


me 2 


so 


dh 

(1tl 


=  1  + 


1 

fL 

c 2 

fN 

- 

(9) 


If  the  receiver  is  a  ground  fixed  station,  then  the  underlined  term  is  constant  for  a  given  station 
and  can  be  handled  separately.  Nevertheless,  we  will  carry  the  term  along. 

For  an  elliptical  orbit 


rT  =  aT  (\  -  cT  cos  Et(tr  (10) 

where  Et(tr  )  is  the  eccentric  anomaly,  er  is  the  eccentricity  of  the  satellite  orbit,  and  aT  is  the  semi- 
major  axis.  Et  is  written  as  a  function  of  tr  since  the  ground  station  clock  is  used  as  a  standard. 
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When  we  solve  Equation  10  for  —  and  substitute  into  Equation  9 

rT 


1 

M  .  R 

c2 

r  1 

_R 

2 aTc2 

+  eT  cos  Et(jr) 
-  eT  cos  Et(t  ) 


but  since 


3  +  e  cos  E  4  e  cos  E 

- 3  + - 

1  -  e  cos  E  1  -  e  cos  E 


Equation  1 1  can  be  rewritten  as 


rr  1  n  3n  vr 

—  =  1  +  —  — - —  + — 


up  2MeJ  cos  Et(tb  ) 


rR  2aT  2  aRc2  1  -  eT  cos  Et(tr  ) 


The  problem  now  is  to  find  how  far  out  of  synchronization  the  clocks  become  after  a  long  time 
interval  (of  the  order  of  seconds,  minutes,  or  even  hours).  To  do  this,  we  integrate  Equation  12 

1  [”m  3m  vi]  2fxeT  fees  Et(tr  )  drR 

tt  tr  c2  r ^  2 aT  2  Tr  <*tc2  J 1  -  eT  cos  Et(tr  )  K 

where  K  is  a  constant  of  integration. 

Kepler’s  equation  is 


where 


n(rR  '  T*  0 )  =  Et^tr  >  "  eTsin  EAtr  > 


n  =  M1/2/ar3/2 


a  3 /2 

^ tr  ^1 12  (*  "  ercos  ET(TR)JdE '7 


Ll  •  I  V  s.  t _ V  .  *■  .  *.  -  n..  \  »  t  >■  t1-  A 
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If  we  substitute  Equation  14  into  Equation  13  and  carry  out  the  integration,  we  get 

1  It  3(1  "l  aTV‘eT  .  ..  ,  ,  ,  „ 

Tr‘rR  *  ~T  +“  rR - 75 - sm£r(r„)  +  K. 


If  desired,  we  can  evaluate  K  in  terms  of  the  time  when  the  clocks  were  synchronized. 
If  tt  =  tr  =  t0,  then  from  Equation  15 


2/xw  a^ej. 


sin  ET(rn) 


1 

V  3/i  (  vr 

c2 

rR  2_ 

For  the  time  of  emission  of  a  pulse,  Equation  1 5  gives 


1 

At 

3At _ 

u! 

c2 

Sr 

2 aT 

2 

2 -  sin  ET(rRe)  +  K 


where  rTg  is  the  time  of  emission  by  the  satellite  clock  and  tr  e  is  the  time  of  emission  by  the  station 
clock.  If  we  let  rRr  be  the  reception  time  by  the  station  clock,  the  uncorrected  range  is  obtained  by 
using  the  satellite  clock  reading  TJe  for  the  emission  time 


f>u=C^R  ~TT  > 


and  the  corrected  range  is 


P=C(X"  T*e) 


where  both  times  are  read  on  the  station  clock. 

From  Equation  17,  the  corrected  range  then  is 

1  n  3#i  VR 

p  =  p  +___  +_  T  . 
“  c  rR  2 aT  2  n  ■ 


ft  ^  p 

-M 


sin  Et(tr  )  +  cK 


where  cK  is  a  range  bias  due  to  relativistic  effects. 
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If  the  user  has  an  orbit  determination  program  that  uses  range  difference  data,  he  can  write 
Equation  18  in  the  form* 


Ap  =  Apu 


2p* 


c 


sin  (Et  +  A Et)  -  sin  ET 


(19) 


where  Er  is  the  eccentric  anomaly  at  the  time  of  emission  of  the  first  pulse  according  to  the  station 
clock. 

Also,  if  the  user  has  already  integrated  an  orbit  such  that  he  has  satellite  position  r  and  velocity  v 
at  various  times,  he  can  avoid  the  calculation  of  ET  in  Equation  18  and  write  the  periodic  relativistic 
correction  as 


2p'/i  arVleT 


sin  E. 


(20) 


This  relation  was  called  to  the  author’s  attention  by  John  T.  Carr  of  the  Space  and  Surface  Systems 
Division,  Strategic  Systems  Department,  NSWC,  Dahlgren,  Va.  The  derivation  of  Equation  20  is  given 
in  Appendix  A. 

Note  that  the  underscored  term  in  all  equations  throughout  is  constant  and  can  be  handled  in 
various  ways.  For  example,  since  the  underscored  quantity  in  Equation  17  causes  the  GPS  clock  to 
run  fast  compared  with  a  station  clock,  prior  to  launch  the  satellite  clock  is  purposely  set  low  to 
10.22999999545  Mhz  or  4.45  x  ICT10  low  in  frequency  relative  to  a  nominal  10.23  Mhz.s  This  in 
effect  removes  the  underscored  term  and  leaves  only  the  periodic  relativistic  term. 


SATELLITE  TO  SATELLITE 

For  the  case  of  satelhte-to-satellite  position  determination,  e.g.,  a  host  vehicle  using  GPS  satellites 
to  navigate  via  GPSPAC,  we  write  Equation  15  for  each  satellite  using  subscript  G  for  GPS  and  H 
for  the  host  vehicle.  We  still  use  R  for  the  ground  station  whose  dock  serves  as  a  standard. 

From  Equation  I  5  then 


Tt;  e  =rRe  + 


p_  3p  ” R 

r  2d  2 


2uv‘  i  .  e 
G  G 


sin  E(.(tr  )  +  A." 


tii  tr  + 

r  r 


1 

P 

3p 

2f'4  V’  ell 

r2 

Ir 

la„ 

T*r 

sin  Eh(tr  )  +  Kf/ 


(21) 


(22) 


where  the  subscripts  e  and  r  mean  emission  and  reception,  respectively,  with  K(.  and  A;/  corresponding 
to  K  of  Equation  1  5. 

*In  going  from  Equation  18  to  Equation  llt,  we  have  used  =  TRr  -  — ,  Ar#  ,  =  A and  dropped 
terms  of  order  1  jc2 .  <’  <’ 


1 
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When  we  use 


1  +  ~2 

C"4 


and 


M  3/a  +^R 


r„  la 


R  (, 


M_  3m  ^r_ 


r„  la „ 


l+? 


JL  llL+Ji 


r„  la 


R  ““A/ 


M  3/a 


rD  la 


R  II 


Equations  21  and  22  can  be  inverted  to  Ol  —  1  to  give 


(?) 


X  TCf  "  ,.2 


M  3m 


y  2 

R  G 


2mw  acwcG 

TCe  +  c2  sin  EG^TrJ  ~ 


T«r  =  T"r 


1 

M 

3m 

n  ^  n  ^  P 

-*1  UH  LH 

c1 

Sr 

‘ 

d 

X  +  c2 

The  uncorrected  range  is 


(21') 


sin  £h(tr  )  -  Kh.  (22') 


Pu  =  c<th. 


tg  > 
<? 


and  the  corrected  range  is 


p  =  c<x "  ■'V 

where  tr^  is  the  time  the  pulse  is  emitted  by  the  GPS  satellite  and  tr  r  is  the  time  of  reception  by  the 
host  vehicle,  each  measured  on  the  earth  station  clock. 


If  Equation  21'  is  subtracted  from  Equation  22'  and  multiplied  by  c 


1 

M 

3m 

i 

M 

3m 

1 

"a* 

• 

P  '  Pu  +  t. 

Sr 

X-  +  -J 

ro 

e 

C 

2m* 


a.”  eH  sin  £„<r„  >-  <V;  sin  lir{r(.  ) 


'll 


-  c  (Kff  -  Ka  ). 


(23) 
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1 

In  Equation  23,  the  station  clock  has  been  eliminated  by  using  tr  =  th  +  O  (— )  and  tr  =  r(; 

j  r  r  C  e 

+  O  (— )  from  Equations  21  and  22. 
c 1 

Equation  23  was  developed  keeping  both  of  the  underscored  terms  in  Equations  21  and  22.  If 
the  GPS  clock  frequency  was  set  low  before  launch  as  discussed  at  the  end  of  the  previous  section, 
then  in  effect  the  first  underscored  term  in  Equation  23  is  removed  and  we  have 


p  =  pu¬ 


la 

la 


H 


VR 
+  — 

1 


2p' 


'//. 


| aH  h  CH  S'n  ^ llT  1!  ^  ~aG  Vl  eG  S*n  ^ G  (TG  g  C^H  ~KG  ^ 


(24) 


leaving  one  cumulative  time  effect.  Also  KG 1  ^  KG  due  to  vanishing  of  the  underscored  term  in 
Equation  16. 


RELATIVISTIC  CORRECTION  TO  DOPPLER  DETERMINATION  OF  RANGE  DIFFERENCE 


In  Appendix  C,  Equation  C-6  gives  the  general  relativistic  Doppler  effect  to  be 


1  + 


2  <P{rT)  v2r 


1  + 


v\ 


[ i  +  $  •  — 

\  c  ) 


(25) 


and  Equation  C-5  gives 

p 

1  -  ~  (26) 

c 

where  subscript  T  refers  to  transmitter  (any  satellite),  subscript  R  refers  to  receiver  (earth  fixed 
dp 

station),  and  p  =  — 
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The  ratio  of  radicals  is  the  same  as  in  the  Relativistic  Correction  to  Range  section,  so  by  the  same 
procedure  used  to  obtain  Equation  12,  Equations  25  and  26  above  become 


1  fa  3  u  vl\  2neT  (  cos  Et(tr) 


rR  2aT  2  )  aTcl  \1  -  eT  cos  l t(tr  ) 


1 


If  we  multiply  the  terms  in  Equation  27  together  and  neglect  terms  of  higher  order  than  1/c2,  this 
becomes 


fr  / n  3p  !  ^T^eT  (  COS  Et(tr)  \  fr  dp 

c2  \rR  2 aT  2  J  aTc 2  \1  -  cos  Et(tr  )f  c  drR  ' 


fR  fT  + 


The  underscored  term  in  Equation  28  can  be  seen  to  be  a  constant  frequency  bias  that  could  be 
solved  in  the  computer  orbit  determination  program.  In  fact,  this  procedure  has  been  followed  over  the 
years  in  Doppler  tracking  of  satellites;  i.e.,  this  relativistic  frequency  bias  has  been  absorbed  into  any 
other  frequency  biases  caused  by  other  effects.  By  handling  the  constant  relativistic  effect  in  this 
manner,  the  only  remaining  relativistic  effect  is  the  periodic  term  caused  by  the  nonzero  eccentricity. 

In  the  following,  we  will  carry  the  constant  term  along  but  continue  to  underscore  it  for  identifica- 


In  tracking  the  satellite,  we  have  a  station  oscillator  of  frequency,  fs ,  with  the  frequency  chosen 
so  that  fs  >  fRmax.  Then,  if  we  mix  fs  with  the  incoming  signal  fR ,  the  beat  frequency  is 


dN 

drR  ~  fs 


so  when  Equation  28  is  subtracted  from  /„ 


dN  _  /  / 3^i  +^V  ^eT  /  cos  ^t^tr  ^  \  |  [t  dp 

dTR  °2  \R  20 T  2  /  aTc2  V1  "  eT  COS  EATR  ^  I  C  dtR 


where  A f  ~  fs~  fT- 


•  -• 
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If  we  integrate  both  sides  of  Equation  29  from  tr  to  tr  +  A tr  and  let  N  be  the  number  of 
Doppler  counts  taken  over  the  time  interval  Ar^ ,  we  get 


N 


T-  Ai*  ^ T  /m  3p  Vr\  a  T^eT  f 

TR  °2  I*  ^  +  2 '  TR  +  °T^  -I  r 


tr  + 
cos 


A  r 


Ej\tr  ) 


1  -  eT  cos  Et(tr  ) 
R 


dTr 


T  A 

+  —  A  p. 

c 


(30) 


As  in  Equation  13,  the  integral  above  is 


2  fTnVlaTv‘eT 

- — -  [sin  ( Et  +  A Et)  -  sin  ET  ]. 

c 

so  if  Equation  30  is  multiplied  by  —  and  solved  for  A p 

fT 

c  r  *|  1  (  n  3u  uj?\  2 nv,aJ.VleT 

^=-fv-A/Ar„]  +-  ^  +  -  — ; -  [sin  (£r  t  A£r)  -  sin  Ar]  .  (31) 


In  Equation  31,  the  first  term  is  the  usual  nonrelativistic  Doppler  result  corresponding  to  Apu 
(uncorrected  range  difference),  the  second  term  is  the  constant  portion  of  the  relativistic  effect,  and  the 
third  term  is  the  periodic  relativistic  effect.  Thus,  Equation  31  above  for  Doppler  data  corresponds 
to  Equation  1 9,  which  was  obtained  by  differencing  two  range  measurements. 
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APPENDIX  A 

AN  ALTERNATE  FORM  FOR  THE  PERIODIC  RELATIVISTIC  EFFECT 


From  Equation  18,  the  periodic  relativistic  effect  (PRE)  is 

fro 


PRE  = 


2vfrcrrie  ... 


sin  E. 


(A-l) 


The  elliptic  orbit  equation  is 


and  Kepler’s  equation  is 


r  =  a  (1  -e  cos  E) 


(t- t0)=*  E-eem  E. 

Differentiating  Equations  A-2  and  A-3  with  respect  to  time 

r  =  ae  sin  E  E 

— =  (1  -  e  cos  E)  E. 

3/2  '  ' 


(A-2) 


(A-3) 


(A-4) 

(A-5) 


When  Equation  A-5  is  solved  for  E,  substituted  into  Equation  A-4,  and  Equation  A-2  is  used  to  eliminate 
the  quantity  ( 1  -  e  cos  E),  one  obtains 


rr  =  n^gfre  sin  E. 


Since  r2  =r  •  r,  then  r  f  =  r  •  r  =  r  •  jy,  so 


PRE  =^-~ ^  sin  E  =  ~r 


(A-6) 


(A-7) 
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APPENDIX  B 


THE  CLASSICAL  DOPPLER  EFFECT 


FIGURE  B-l.  RECEIVER  AT  REST,  TRANSMITTER  MOVING  AWAY 

When  T  is  at  position  xQ,  we  imagine  an  event  occurring  at  its  origin  (e.g.,  a  light  flash)  at  time  t\  by 
T'b  clock.  Since  the  light  travels  at  velocity  c,  when  R  sees  the  event  its  clock  reads 

f,  =  t{  +-^  (B-l) 

a  later  time  because  of  the  transmission  time  xQ/c  (Figure  B-l). 

At  a  later  time,  t2,  another  event  takes  place;  R  sees  this  event  when  its  clock  reads 

t +  (B-2) 

c  c 

The  additional  distance  of  travel  is  vf  (t2  -  t{),  not  vT  (t2  -  f, )  since  the  former  is  the  true  distance 
that  T  has  traveled  in  the  actual  time  interval  between  the  events,  while  t2  -  tt  is  the  apparent  time 
interval  according  to  R. 

When  Equation  B-l  is  subtracted  from  Equation  B-2 

fj-*i  (ti  -  t,') 


c  it:  ptfWTn,  n 


(WWWW 


NSWC  TR  83-55 


or  let  A  tR  =  t2  -  ft ,  Af  T  =  t2  -  t2  and  we  have 


AtR  =  Atr  (1  +— )•  (B-3) 

c 

This  is  the  difference  in  time  intervals  because  of  the  motion  changing  the  transmission  times.  The 
clocks  are  assumed  to  be  identical  and  running  at  the  same  rates;  i.e.,  no  relativistic  effects. 


Equation  B-3  is  (classically)  true  for  any  pair  of  events;  they  do  not  have  to  be  periodic.  If  the  events 

are  periodic  (such  as  the  vibrations  of  an  electron  producing  a  light  wave),  then  ^  _ _ !_  ,  -■  _  1 

so  using  Equation  B-3  R  T  At T 


f  -  T 

f*  "  1  +Vt  <B‘4) 


FIGURE  B-2.  TRANSMITTER  AT  REST,  RECEIVER  MOVING  AWAY 


In  the  same  manner  as  before  (Figure  B-2) 

(B-5) 

and 

<»“  *2  +—+-!£-(** -ft').  (B-6) 

Note  that  one  might  be  tempted  to  write  for  the  additional  transmission  distance  in  Equation  B-6 
vK  (t2  -  1 1 )  since  t2  ~tx  is  the  true  time  interval  between  the  events.  R,  however,  which  is  in  motion, 
has  been  traveling  for  a  time  t2  -  during  the  interval  in  which  it  has  observed  both  events.  The  addi¬ 
tional  distance  is  then  vR  (t2  ~t2). 

If  Equation  B-5  is  subtracted  from  Equation  B-6 

t2  -  ti  =  t2  - 12  +_*L  (tj  -  f/), 
c 


and  put  Atn  =  t2  -  f,',  A tT  =  1 2  -  f, ,  then 
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At/?  =  A*r  +ir  At R 


or 


A  tD  = 


Af, 


c 


(B-7) 


Again,  assuming  perfect  clocks,  this  is  the  difference  in  time  intervals  because  of  the  transmission 


time. 


If  the  events  are  periodic,  fR  =  ,  fT  =  so  from  Equation  B-7 

R  T 


fR  =fTa-vJL). 

c 


(B-8) 


To  summarize  to  this  point 


ft j  -  ~  for  R  at  rest,  T  moving  away, 


1+^ 


(B-4) 


fR  ~  fj  for  T at  rest,  R  moving  away. 


(B-8) 


FIGURE  B-3.  TRANSMITTER  AND  RECEIVER  MOVING  IN  ARBITRARY  DIRECTIONS 


’  vl 


To  be  more  general,  both  Equation  B-4  and  Equation  B-8  assumed  R  and  T  were  separating,  so  in 
those  equations  -vR  -*  p  •  vT ,  Vj-  -*■  p  •  Vf  where  p  is  a  unit  vector  directed  from  R  toward  T  (Figure  B-3). 
So  Equation  B-4  and  Equation  B-8  can  be  written 

fR  =  _  R  at  rest  (B-9) 

1+  P'It 

c 

fR  afT(l+Pm  )  T  at  rest.  (B-10) 
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To  put  both  R  and  T  in  motion,  we  imagine  a  “virtual”  observer  V  at  rest  in  inertial  space  on  the  line 
joining  R  and  T.  Then,  since  the  observer  is  at  rest,  the  frequency  it  receives  from  T  according  to  Equa¬ 
tion  B-9  is 


and  the  frequency  received  by  R  transmitted  by  V  is  by  Equation  B-10 


Since  V  is  an  empty  point  in  space,  the  signal  simply  passes  through  so 


/t  ft} 
ly  Ry 

When  we  substitute  Equation  B-ll  into  Equation  B-12  using  Equation  B-13,  we  have 


1+P 


1  +  p  •  — — 

c 


(B-ll) 


(B-12) 


(B-13) 


(B-14) 


which  is  the  classical  Doppler  equation  when  both  R  and  T  are  in  motion. 

In  Doppler  tracking  of  satellites,  we  write  the  Doppler  equation  as 

/*’(1‘*)/r  <B1B) 

where  p  is  the  time  rate  of  change  of  the  range  between  the  station  and  satellite.  At  first  glance,  it  might 
appear  that  in  Equation  B-14  we  are  neglecting  the  receiver  (station)  velocity  and  writing  p  •  tp  -  p  so 
that  to  first  order  Equation  B-14  is 


This  is  not  the  case! 

To  derive  Equation  B-15,  we  write  the  range  vector  as 

P  -rrtteJ-rR  (*r)  (B-16) 

where  te  -  time  emitted  and  tr  =  time  received.  The  time  emitted  is  obtained  from 

From  Equation  B-16 

A  dFr(te)  dtp  (*r) 
p  dtr  ~  dtr 
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u  =  drT  (te)dte_drK  ( tr )  . 
dte  dty  dtr 

From  Equation  B-17  =  1-—  so  (B-18)  is 

dtr  c 


(B-18) 


(B-19) 


Define  a  unit  vector  p  directed  from  the  station  toward  the  satellite.  Then  if  we  take  the  dot  product  of 
p  with  p\  we  obtain  from  Equation  B-19 


From  Equation  B-20 


P  =  P 


P 


P 


P 


p  •  1jt  -  p  -  1jr 

1  +p  •  Zl 


so 


I 


f 


1+p  •  f 

i  ^ T 

1  +p  . 

c 


(B-20) 


(B-21) 


When  Equation  B-21  is  substituted  into  Equation  B-14,  the  Doppler  equation  is  put  in  the  form  of 
Equation  B-15.  It  is  clear  then  that  Equation  B-15  is  an  exact  expression  of  the  classical  Doppler 
equation.  There  is  no  dropping  of  p2  /c2  and  higher-order  terms,  and  the  station  is  not  considered  at  rest. 
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APPENDIX  C 

THE  RELATIVISTIC  DOPPLER  EFFECT 


The  proper  time  in  any  system  is  given  by  Equation  5.  We  can  then  write  for  transmitter  and 
receiver 


dr 


T 


drD 


(C-l) 


(C-2) 


Note  that  here,  unlike  Equations  6  and  7,  we  have  put  subscripts  on  the  coordinate  times.  The  reason 
is  that  in  the  earlier  section  of  the  report  our  goal  was  to  compare  clock  rates  assuming  that  the  trans¬ 
mission  time  was  taken  into  account. 


The  Doppler  effect,  however,  arises  because  of  the  transmission  time,  as  was  shown  in  Appendix  B. 
Toward  our  goal  of  deriving  the  relativistic  Doppler  effect  then,  tT  in  Equation  C-l  is  the  coordinate  time 
when  the  satellite  emitted  the  signal,  and  tR  in  Equation  C-2  is  the  coordinate  time  when  it  is  received  by 
the  station. 


The  proper  frequency  is  fR  = 


1 

drR 


»  f t 


—  so  if  we  divide  Equation  C-l  by  Equation  C-2 


Lb. 

fT 


V,  .  2<Mrr)  t>r» 

1  c*  ~c* 

V 1  +  “pli? 

C 1  C 2 


(C-3) 


In  order  to  obtain  an  expression  for  —  £. ,  we  write  from  Equation  C-2 

dtR 


r 


NSWC  TR  83-55 


so 


dtT  drRe 
dtR  drRr 


Since  the  corrected  range  is 


P  =  c(TRr~  TRe ) 


then 

r  —  t  —  — 
Re  Rr  c 

so 


dtT  drR e 


dtR  dTRr 


From  Appendix  B,  Equation  B-21,  we  see  that  Equation  C-4  can  be  written  as 


(C-4) 


(C-5) 


so  Equation  C-3  becomes 


UL  J 
ff 


^JVi  .  2*(rr>  _^/L 
c2  c2 


(C-6) 


Equation  C-6  is  the  general  relativistic  Doppler  effect.  If  i?  and  T  are  in  free  space  far  from  any 
source  of  gravitation,  4>(/y)  =i<b(rji )  =  o  and  Equation  C-6  reduces  to  the  special  relativistic  Doppler 

effect.  If  also  <1  and  <1,  then  Equation  C-6  reduces  to  the  classical  Doppler  effect. 
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